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occupancies for static properties of nuclei. 
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' Abstract. The possibility to use functionals of occupation numbers and natural orbitals for inter- 
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INTRODUCTION 

The nuclear many-body problem of N interacting nucleons can be solved exactly only 
in very specific cases or for very small particle numbers. This is due to the large number 
of degrees of freedom involved in such a complex system. Let us for instance consider 
particles interacting through n Hamiltonian written as 

; H = Y i t ij a i a j + J Y,VijU<4'tfaia k +-~ (1) 
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^! . Then the exact ground state energy can be written as 
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where fj-' = (ataj), y^L = (af a- aia^), ... denote the one-, two-, ... body density 
matrices that contain all the information on the one-, two-... body degrees of freedom 
respectively. A natural way to reduce the complexity of this problem is to assume that 
at a given level, the k— body (and higher-order) density matrices becomes a functional 
of the lower-order ones. This is what is done for instance in the Hartree-Fock (HF) 
approximation where all fc-body density matrices (with k > 2) become a functional of 
y^ 1 ). Unfortunately, the HF theory applied to the nuclear many -body problem in terms of 
the vacuum Hamiltonian is a poor approximation and Many-Body theories beyond HF 
are necessary. 

The introduction of Energy Density Functional (EDF) approaches in the 70's was a 
major breakthrough (see for instance [1] for a recent review). In its simplest form, the 
EDF formalism starts with an energy postulated as a functional of the latter being 
built out of a Slater Determinant. Then the ground state energy is obtained by minimizing 



the energy with respect to y( l \ i.e. 
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Parameters are generally adjusted on specific experimental observations and therefore 
encompass directly many-body correlations. Current EDF uses a generalization of eq. 
© obtained by considering quasi-particle vacua as trial states. By making explicit use 
of symmetry breaking, such a functional called hereafter Single-Reference (SR-) EDF 
is able to account for static correlation associated with pairing and deformation. Actual 
SR-EDF takes the formQ: 

^Exact 

where K denotes the anomalous density. To restore symmetries and/or incorporate dy- 
namical correlations, guided by the Generator Coordinate Method (GCM), a second 
level of EDF implementation, namely Multi-Reference (MR-) EDF is introduced. Re- 
cently, difficulties with the formulation and implementation of have been encountered in 
MR- EDF. A minimal solution has been proposed in ref. J^HHH]. Besides these problems, 
the authors of ref. $2^ have pointed out the absence of a rigorous theoretical framework 
for the MR EDF approach. At the heart of the problem is the possibility to break sym- 
metries in functional theories and then restore them using configuration mixing. This 
issue needs to be thoroughly addressed in the future. 

In this context, it is interesting to see if extensions of the functional used at the SR- 
EDF level can grasp part of the effects that for standard functionals require the MR level. 
It is worth realizing that, in the canonical basis for which = \ (pj)rii((pj\, we have 

= <f C or[{<Pf,"/}] = \ £ VffJjV^O- " "/) \/ ^ ~ *j) , (5) 
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and therefore, the energy can be regarded as a functional of natural orbitals <p,- and oc- 
cupation numbers As a matter of fact, for electronic systems, Gilbert has generalized 
the Kohn-Sham theory and shown that the exact energy of a system can be obtained 
by minimizing such a functional [5] leading to the so-called Density Matrix Functional 
Theory (DMFT). The possibility to consider occupation numbers as building blocks of 
the nuclear energy functional has recently been discussed in ref. [0, @]. Two levels of 
theory can be developed along the line of Gilbert's idea (i) either, functionals in the 
strict Gilbert framework can be designed. In that case, since the density identify with the 
exact density at the minimum, it should respect all symmetries of the bare Hamiltonian. 
(ii) or we exploit the concept of symmetry breaking. In the latter case, similarly to the 
SR-EDF, strictly speaking we cannot anymore rely on the theorem, but we may gain 
better physical insight with relatively simple functionals. 



Note that the denomination "mean-field" or the separation into a "mean-field" like and "correlation" like 
is completely arbitrary since, as we mention previously, the so-called "mean-field" part already contains 
correlation much beyond a pure Hartree-Fock approach. 



APPLICATION TO THE LIPKIN MODEL AND DISCUSSION 



The descriptive power of DMFT is illustrated here in the two-level Lipkin model HlOll . 
In this model, the Hartree-Fock (HF) theory fails to reproduce the ground state en- 
ergy whereas configuration mixing like Generator Coordinate Method (GCM) pro- 
vides a suitable tool da |9Q. Therefore, the two-level Lipkin model is perfectly suited 
both to illustrate that DMFT could be a valuable tool and to provide an example of 
a functional for system with a "shape" like phase-transition. In this model, one con- 
siders N particles distributed in two N-fold degenerated shells separated by an en- 
ergy e. The associated Hamiltonian is given by H — £Jq — j(J + J + + /_/_) where V 
denotes the interaction strength while Jq, J± are the quasi-spin operators defined as 

Jo = 3 Lp=i (ci, p c+,p ~ c -,p c -, P )> J + = Ip=i c +,p c -, P and J - = 4- C +, P and C -, P are 
creation operators associated with the upper and lower levels respectively. Due to the 

specific form of the Lipkin Hamiltonian , simply writes in the natural basis as 
7 (1) = Lp=\ {l<Po,/?W(<Po, P | + |<Pi,p)«i(<Pi, P |} with ni = (1 -n ). Introducing the an- 
gle a between the state \—,p) and \(po,p), leads to the following mean-field functional 
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4*f({<P*,W) =£M<*,no) = -|^v{cos(2a)(2n -l) + |sin 2 (2a)(2no-l) 2 }. (6) 

where X = V(N — l)/£- This expression is easily obtained by generalizing the Hartree- 
Fock case (recovered here if hq = 1). The main challenge of the method is to obtain an 
accurate expression for <fcor- To get the functional, clearly identified cases from which 
properties of the functional could be inferred have been used rflin . namely the N = 2 case 
and the large N limit. In the two-particles case, the correlation energy can be analytically 
obtained and reads 

<fcor 2 (cMo) = -2V| sin 2 (2a)n (l -no) + (sin 4 (a) +cos 4 (a)) a/«o(1 -»o)}(7) 

A simple extension of the N = 2 case for larger number of particles is to assume that each 

pair contributes independently from the others leading to Sq ot = [N(N — \ )/2]$ Cm . 
However, such a simple assumption leads to a wrong scaling behavior in the large 
N limit. Indeed, in this case, <^ r °^ N 2 as N tends to infinity while a N 4 ^ 3 scaling 
is expected H12I1 . To obtain the correct limit, a semi-empirical factor 7}(N) can be 
introduced such that 

N>3 N(N — I) N=2 

S Cor (a,n ) = t?(A0 - V Cor (a, no), (8) 

with T}(N) = cN~ 2 / 3 . The value c = 1.5 has been retained using a fitting procedure. 
Examples of results obtained by minimizing the functional given by Eqs. © and ([8]) are 
shown in Fig. [TJ for different particle numbers and interaction strengths. In all cases, a 
very good agreement, much better than the HF case is found. 

The Lipkin example suggests that DMFT can be a valuable tool for describing ground 
state of a many-body system when symmetry breaking plays a significant role. The 
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FIGURE 1. Exact ground state energy (solid lines) displayed as a function of % for N — 5 to 20 resp. 
from top to bottom. In each case, the corresponding HF (dashed line) and DMFT (filled circle) minimum 
energy are shown. The DMFT calculation is performed using the mean-field and correlation energy resp. 
given by Eq. © and Eq. © with rj (N) = 1 .5 2/3 (Adapted from ifTl l. 

functional designed here is exact only in the N = 2. Note that the functional proposed 
here breaks signature symmetry and therefore enters into the level (ii) of functional 
discussed in the introduction. The Lipkin model is however rather schematic and cannot 
be used as a guidance for realistic situations. The possibility to design a new accurate 
functional for nuclei remains a challenging problem. 
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